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Abstract 

Effects of correlated disorder on wave localization have attracted considerable interest. Motivated 
by the importance of studies of quantum transport in rough nanowires, here we examine how 
colored surface roughness impacts the conductance of two-dimensional quantum waveguides, using 
direct scattering calculations based on the reaction matrix approach. The computational results 
are analyzed in connection with a theoretical relation between the localization length and the 
structure factor of correlated disorder. We also examine and discuss several cases that have not 
been treated theoretically or are beyond the validity regime of available theories. Results indicate 
that conductance properties of quantum wires are controllable via colored surface disorder. 
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I. INTRODUCTION 

Ever since Anderson's model of electron transport in disordered crystals [lj, wave lo- 
calization in disordered media has attracted great interests due to their universality. For 
example, two recent experiments directly observed matter-wave localization in disordered 
optical potentials using Bose-Einstein condensates One of the most known results 

from Anderson's model is that in one-dimensional (ID) disordered systems, the electron 
wavefunction is always exponentially localized and hence does not contribute to conduc- 
tance for any given strength of disorder. Note however, this seminal result is based on the 
strong assumption that the disorder is of the white noise type. If the disorder is colored due 
to long-range correlations, then a mobility edge may occur in one-dimensional systems as 

weii y. 

Quantum transport in nanowires is of great interest due to their potential applications 
in nanotechnology. In addition to the possibility of ballistic electron transport, quan- 
tum nanowires are found to show many other important properties. In particular, silicon 
nanowires can have better electronic response time [5fl as well as desirable thermoelectric 

n 

properties [6J. It is hence important to ask how the nature of surface disorder of quantum 
wires, modeled by quantum waveguides in this study, affects their conductance properties. 

Remarkably, if the surface scattering contribution is weak, then it is possible to map 
the conduction problem of a long two-dimensional (2D) rough waveguide to that of a ID 
Anderson model of localization, with the disorder potential determined by the surface rough- 

on 

ness p, LSf] . Initially this mapping was established for one- mode scattering but later it was 
generalized for any number of modes in the transverse direction |9|]. As such, a quan- 
tum wire with white-noise surface disorder will have zero conductance if the localization 
length is much smaller than the wire length. However, in reality the surface disorder of a 
rough quantum wire always contains correlations. As a result it becomes interesting and 
necessary to understand the conductance properties in rough quantum wires with their sur- 
face disorder modeled by colored noise. This has motivated several pioneering theoretical 



studies 
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Under certain approximations the theoretical studies predicted 



localization-delocalization transitions of electrons in 2D waveguides with colored surface dis- 
order. Some theoretical details were tested by examining the eigenstates of a closed system 
with rough boundaries [12]. Moreover, the predicted mobility edge due to colored disorder 
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was recently confirmed in a microwave experiment [13]. 

Using a reaction matrix formalism for direct scattering calculations, here we computation- 
ally study the conductance properties of rough quantum wires with colored surface disorder. 
The motivation is threefold. First, though the dependence of the localization length upon 
the correlation function of surface roughness is now available from theory, how the more 
measurable quantity, namely, the conductance of the waveguide, depends on colored sur- 
face roughness has not been directly examined. This issue can be quite complicated when 
the localization length becomes comparable to the waveguide length. Second, computation- 
ally speaking it is possible to consider any kind of colored surface disorder, thus realizing 
interesting circumstances that are not readily testable in today's experiments. Indeed, in 
our computational study we can create rather arbitrary structures in the surface disorder 
correlation function and then examine the associated conductance properties. Third, di- 
rect computational studies allow us to predict some interesting conductance properties that 
lave not been treated theoretically or go beyond the validity regime of available theories 
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14j • For example, we shall study the conductance properties for very strong sur- 



face roughness, for rough bended waveguides, and for scattering energies that are close to 
a shifted threshold value for transmission. The long-term goal of our computational efforts 
would be to explore the usefulness of colored surface disorder in controlling the conductance 
properties. 

This paper is organized as follows. In Sec. II we describe the scattering model of a 
quantum wire with colored surface disorder. Therein we shall also briefly introduce the 
methodology we adopt for the scattering calculations. In Sec. Ill we present detailed 
conductance results in a variety of one-mode scattering cases and discuss these results in 
connection with theory. Concluding remarks are made in Sec. IV. 



II. QUANTUM SCATTERING IN WAVEGUIDES WITH COLORED SURFACE 
DISORDER 

A. Modeling waveguides with colored surface disorder 

We treat quantum wires as a long 2D waveguide as illustrated in Fig. 1(b). The scattering 
coordinate is denoted x and the transverse coordinate is denoted y. The width of the 
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FIG. 1: Schematic plot of a 2D rough waveguide that models rough quantum wires, (a) The 
generation of a rough surface is illustrated using M = 4 random shifts in the transverse direction, 
(b) One waveguide geometry with a straight upper boundary y = P{x) = 1 and a rough lower 
boundary y = Q(x). Scattering occurs in region I (gray area) and region II denotes the left and 
right leads. Arrows indicate the direction of incoming and outgoing electron waves. A n , B n , C n 
and D n are quantum amplitudes (see Eq. (8)). 

waveguide is denoted w and the length is denoted L. In all the calculations L = lOOw and 
w is set to be unity. That is, we scale all length by the waveguide width. The upper and 
lower boundaries of the waveguide are described by y = P(x) and y = Q(x). The case in 
Fig. 1(b) represents a situation where the upper boundary is a straight line (P(x) = 1) and 
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the lower boundary is rough. As in our other studies of rough waveguides 
a rough waveguide boundary in three steps. First, we divide a rectangular waveguide into 
M pieces of equal length L/M. Second, the end of each piece is shifted in y randomly, with 
the random y-displacement, denoted rj, satisfying a Gaussian distribution. Third, we use 
spline interpolation to combine those sharp edges to generate a smooth curve r](x) for either 
the upper or the lower waveguide boundary. For the sake of clarity, Fig. 1(a) depicts this 
procedure with the number of random shifts being as small as M = 4. In all our calculations 
below we set M = 100. In Fig. 2(a) we show one realization of the surface roughness function 
rj(x). 

The function i](x) may be characterized by its ensemble- averaged mean rj and its self- 
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correlation function C v (x — x'), i.e., 

rj = (rj(x)) = 0, 

(rj(x)T)(x')) = a 2 C v (x - x'), (1) 

where a is the variance of r](x). In the limit of white noise roughness, C v (x — x') is propor- 
tional to 5(x — x'). But more typically, C n (x — x') decays at a characteristic length scale, 
called the correlation radius R c . For our case here, because the randomness is introduced 
after dividing the waveguide into M = 100 pieces, the correlation length R c of the rj(x) we 
construct here is of the order of L/M ~ w. This length scale is comparable to the wavelength 
of the scattering electrons in the one-mode regime. 

One tends to characterize the strength of the surface roughness by the variance a defined 
above. However, in practice it is better to use the maximal absolute value of rj(x), denoted 
l^maxl; to characterize the roughness strength. This is because for strong roughness with 
a given variance, there is a possibility that some of the random displacements become too 
large such that the waveguide may be completely blocked. Recognizing this issue, we first 
set a value of |?7 max | and then, after having generated a roughness function r](x) based on 
spline interpolation, rescale r](x) such that its maximal absolute value is given by \r) max \ . 

The roughness function rj(x) obtained above does not have any peculiar features. There 
are a number of ways to introduce some structure to the correlation function C n {x — x'). 



In Ref. 17| a filtering function method was proposed to produce a power-law decay of 
C v (x — x') from white noise. Here we adopt the approach used in Ref. jis| . which is 
based on the convolution theorem of Fourier transformations. In particular, the discrete 
form of autocorrelation function of rj(x) is defined as 

yr ,mL N Jpy 1 f(n + m)L\ fnL 



n=l 



where m = —N +!,•••, —1, 0, 1, • • • , N — 1, iV is the total number of grid points along 



x, and c is a normalization constant such that C v (0) = 1 19j. In Fig. 2(b) we show 
the autocorrelation function for the surface roughness function depicted in Fig. 2(a). The 
autocorrelation drops from its peak value to near zero at a scale of R c ~ 0.7w, which is 
much smaller than the waveguide length. 

As will be made clear in what follows, it is important to consider the Fourier transform 
of C v (x), i.e., the autocorrelation function in the Fourier space. This important quantity is 




FIG. 2: (a) One realization of the surface roughness function rj(x), with the method described 
in detail in the text, (b) The associated autocorrelation function C v (x). (c) Surface structure 
factor Xri(k) obtained from the C n {x) shown in panel (b). (d) A function p(x) that will be used to 
introduce additional correlations via convolution, (e) The Fourier transform of the p(x) shown in 
panel (d). (f) The structure factor Xf)(k) obtained from a convolution between p{x) and r](x). (g) 
The new surface roughness function fj(x), with correlations that are absent in rj(x). 

denoted x^ik), where k is the wavevector conjugate to x. Using the Fast Fourier transform 
of C v (x), Xv(k) can be evaluated as follows: 



2N 



3=1 



-i2ir(j-l)(m-l) 
2N 



(3) 



where the value of k on the left side is determined by the value of m on the right side via 
k = (2(m/2N) - \){2-kN/2L). In our calculations we choose TV = 1024. Note that x v ( k ) 
is a real function due to the evenness of C v (x). The real function x-qik) is called below 
the structure factor of the surface roughness. Figure 2(c) shows the structure factor Xr,{k) 
obtained from the correlation function shown in Fig. 2(b). 

Additional correlations in the surface disorder can now be generated by modulating the 
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structure factor Xv(k)- Because the structure factor Xrj(k) for a single realization is equivalent 
to the square of the Fourier transform of rj(x), we may imprint interesting structures onto 
Xr/{k) by convoluting r](x) with some filtering function. Consider then the function p'(x) = 
sin(ax)/aa; with a > 0. Its Fourier transform is a step function of \k\, ^ with a height n/a 
and the step edge located at \k\ = a. Consider then a combination of n such functions, i.e., 

\ a n\ X 

where A n , a r n , a l n are predefined parameters. Then the Fourier transform will be 7i/A n if 
a r n > \k\ > a l n or a r n < \k\ < a l n ; and zero otherwise. If we now consider the following 



mL\ v-^ / nL\ f(m — n)L\ 



roughness function |20| 

/; 1 N I ^ r \N I '' V A 

then according to the convolution theorem, we have 

lx#)| 1/2 Hxp(*)llx#)i 1/2 , (6) 

where Xfj(k) is the structure factor for the new surface roughness function fj(x). As such, 
the structure of Xp(k) is directly imprinted on Xv(k). That is, computationally speaking, 
arbitrary modulation can be imposed on the structure factor by filtering out the unwanted 
components and magnifying other desired structure components. Below we apply this sim- 
ple technique to create different kinds of surface roughness correlation windows and then 
examine the conductance properties. In Fig. 2(d) we show one example of p(x). Its Fourier 
transform amplitude, as shown in Fig. 2(e), displays two windows. As shown in Fig. 2(f), 
this double-window structure is passed to Xfj(k) due to Eq. (jSJ). Finally, in Fig. 2(g) we 
show the surface roughness function fj(x), which obviously contains more correlations than 
the old surface roughness function r)(x) shown in Fig. 2(a). 

B. Reaction matrix and scattering matrix 

Here we briefly describe how we calculate the electron conductance of a rough 2D waveg- 
uide as described above. The Hamiltonian for the quantum transport problem is given 
by 

h 2 ( d 2 d 2 \ 



2m* \dx 2 dy' 
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where m* is the electron effective mass and V c (x, y) represents a hard wall confinement 
potential. That is, V r (x, y ) is zero in Q(x) < y < P(x) and becomes infinite otherwise. 



In our early work [15|, ll6j| we formulated such a waveguide scattering problem in detail 
in terms of the so-called reaction matrix method. In the reaction matrix method we first 
expand scattering state in the scattering region (region I, gray area in Fig. 1(b)) in terms 
of a complete set of basis states. The basis states are obtained by transforming the rough 
waveguide into a rectangular one, with the expense of a transformed Hamiltonian with extra 
surface dependent terms. The solutions in the leads (region II, Fig. 1(b)) are given by 

C n ^-D n t—)^1) (8) 
\fk n x/K J w 



for the left and right leads, respectively. Here n is the index for the modes in the transverse 
direction, and the wavevector k n is given by 



2mE /mr\ 2 

where E is the initial electron energy. The scattering coefficients A n , B n , C n and D n in 
Eq. (|S]) are determined by the scattering matrix 5", which relates the outgoing states to the 
incoming states. Specifically, 




(10) 



where the submatrices r and r' denote the reflection matrix and t and t' denote the transmis- 
sion matrix. In the case of one- mode scattering (n = 1) considered below, ki will be simply 
denoted as k, with < kw/n < \/3. The S matrix is related to the so-called i?-matrix in 
the reaction matrix method as follows, 

I 2m -iKRK 

I 2m + iKRK- 1 ' 

where m is maximal number of propagating modes and /2m is a 2m x 2m unit matrix, and K 
is 2m x 2m diagonal matrix with diagonal elements determined by the wavevector associated 



with each scattering channel 15|, ll6(]. Once the S matrix is obtained from the R matrix, 



the conductance is calculated by G = GoTrace(tt'), where Go = e 2 /(2h) is the conductance 
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quanta. Note that in our calculations we include about 10 evanescent modes though we focus 
on the energy regime where only one mode in the ^-direction admits propagation along x. 
As to the number of basis states we use in describing the transformed rectangular waveguide, 
we use 1000 basis states for the x degree of freedom and 4 basis states for the y degree of 
freedom. Such a large number of basis states is for a good description of the scattering wave 
function inside the waveguide, and this number should not be confused with the number of 
propagating modes or evanescent. Good convergence is obtained in our calculations. Note 
also that due to the large number of basis states used in the scattering direction, the Fourier 
transform techniques developed in Ref. [15] is especially helpful. 



III. EFFECTS OF COLORED SURFACE DISORDER ON CONDUCTANCE 

With the mapping between the scattering problem in 2D waveguide and ID Anderson's 
model [3, l8j, early theoretical work established that the localization length L\ oc of the 



2D waveguide problem is given by 



loc " (2A;) 2 ' 



where x(2&) is either the structure factor x r? (2A;) or the new structure factor Xfjfik) after 
a convolution procedure. If L loc > L, a transmitting state is expected and if L\ oc « L, 
then the electron can only make an exponentially small contribution to the conductance. 
As such, one expects transmitting states when the structure factor x(2&) is essentially zero; 
and negligible conductance if x(2A;) is significant and if a is not too small. This suggests 
that the conductance properties can be manipulated by realizing different surface roughness 
functions. 

Equation ffl2l is obtained under a weak electron scattering approximation (Born approx- 
imation). As such, the theoretical result of Eq. (fl2l) may not be valid if a is not small 
as compared with w or if the scattering electron is close to the threshold value of channel 
opening. Another assumption in the theory is that L\ oc should be much greater than R c , 
the radius of the surface correlation function C n {x). 
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FIG. 3: (a) Conductance of rough waveguides vs k = k\ (see Eq. (9)). The upper boundary- 
is flat, i.e., P{x) = 1, and the lower boundary is given by Q(x) = r](x), with a surface disorder 
strength characterized by |r? m ax| = 0.2w and the variance a given by a = 0.0779tu, a = 0.0802-u? 
and a = 0.0773tu. Inset on the right shows the structure factor of the surface roughness in one 
single realization, (b) Same as in panel (a), but two cases with different surface structure factor 
obtained from a convolution approach are plotted, using solid and dashed lines. All conductance 
curves here are averaged over three realizations. 

A. Straight Rough Waveguides 

However, in our computational studies we will examine some interesting cases that are 
evidently beyond the validity regime of the theory. For example, the strength of the surface 
disorder may not be small and the scattering energy may be placed in the vicinity of a 
shifted channel opening energy. 

In Fig. 3(a) we show conductance results averaged over three realizations of a rough 
waveguide, with a flat upper boundary P(x) = 1 and a rough lower boundary Q(x) = rj(x). 
The strength of the surface disorder is characterized by |?7max| = 0.2w. Due to our procedure 
in generating a fixed |?7 ma x| = 0.2w, the variance a of the surface function in each single 
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realization of the surface function will change slightly. For the three realizations used in 
Fig. 3(a), a = 0.0779u>, 0.0802w and 0.0773w. As is clear from Fig. 3(a), there exists a 
threshold k ~ 0.6n/w beyond which the system becomes transmitting (this threshold will be 
explained below). In the transmission regime the conductance shows a systematic trend of 
increase as the wavevector k increases. The inset of Fig. 3(a) shows Xr?(2/c), one key term in 
Eq. (fl2|) . The characteristic magnitude of Xr?(2A;) for the shown regime of k is ~ 0.3. Using 
Eq. ffl2l) . one obtains that the localization length L\ oc is comparable to L — lOOw. This 
prediction is hence consistent with our computational results that demonstrate considerable 
transmission. 

Next we exploit the convolution technique described above to form new rough surfaces 
described by fj(x). In particular, the inset of Fig. 3(b) shows two sample cases with dis- 
tinctively different surface structure factors. In one case (dotted line) x^(2fc) has significant 
values in the interval 0.67 < kw/n < 0.8. Indeed, during that regime the value of Xrj(2&) is 
many times larger than the mean value of Xrii^k) in the case of Fig. 3(a). In the other case 
(solid line) Xvi^k) is large only in the regime of 0.75 < kw/n < 0.9. For these regimes, the 
theory predicts that the localization length to be much smaller than the waveguide length 
and hence vanishing conductance. This is indeed what we observe in our computational 
study. As shown in Fig. 3(b), either the dotted or the solid conductance curve display a 
sharp dip in a regime that matches the main profile of x^(2/c). 

In addition, similar to what is observed in Fig. 3(a), Fig. 3(b) also displays a transmission 
threshold. Take the dotted line in Fig. 3(b) as an example. For kw/n < 0.55, there is no 
transmission at all, even though x^(2/c) in that regime is essentially zero. This suggests that 
this threshold behavior is unrelated to surface roughness details. Rather, it can be considered 
as a non-perturbative result that is not captured by Eq. ffT2]) . To qualitatively explain 
the observed threshold, we realize that due to the relatively strong surface roughness, the 
effective width of the waveguide decreases and as a result, the effective mode opening energy 
increases [l5[. For |r7 max | = 0.2w, we estimate that the effective width of the waveguide is 
given by w e s = w — \r] mauX \ = 0.8w. Hence, the corrected mode opening energy E is now 
given by (h 2 /2m*)(ir /0.8w) 2 . Using Eq. ([9]), this estimate gives that, regardless of the 
surface roughness details, the threshold k value for transmission is ~ 0.75ir/w, which is 
close to what is observed in Fig. 3. Such an explanation is further confirmed below. This 
also demonstrates that the maximal value of \r){x) \ is an important quantity to characterize 
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FIG. 4: Conductance in rough waveguides with increasing strength of surface roughness. The 
waveguide geometry is similar to what is considered in Fig. 3, but with a different structure factor. 
\Vmax.\ = 0.01 and a = 0.0046u; in (a), |r/ max | = 0.1 and a = 0.0400k; in (b), |r/ max | = 0.2 and 
a = 0.0779u; in (c), and |r? m ax| = 0.3 and a = 0.1255k; in (d). Inset shows the surface structure 
factor used in (a)-(d). 

the surface roughness strength. Of course, the exact dependence of the effective waveguide 
width upon r] max is beyond the scope of this work [211 ]. 

The results in Fig. 3 show that even when the surface roughness is strong enough to 
significantly shift the threshold energy for transmission, the surface structure factor may 
still be well imprinted on the conductance curve. Moreover, the resultant windows of the 
conductance curves in Fig. 3 are seen to match the location of the structure factor peak. 
Nevertheless, one wonders how such an agreement might change if we tune the strength of 
the surface roughness. 

To that end we examine in Fig. 4 four scattering cases with increasing roughness strength, 
with |?7 m ax| = 0.01k; and a = 0.0046k; in Fig. 4(a) (representing a case with quite weak 
surface roughness), |r/ max | = 0.1k; and a = 0.0400k; in Fig. 4(b), |r/ max | = 0.2k; and a = 
0.0779k; in Fig. 4(c), and |r7 max | = 0.3k; and a = 0.1255k> in Fig. 4(d) (representing a 
case with very strong surface roughness). The main profile of the structure factor is also 
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shifted closer to the threshold regime observed in Fig. 3. For the case with |r] max | = 0.3w, 
the theory based on the weak roughness assumption is not expected to hold. Indeed, In 
Fig. 4(d) the transmission threshold is right shifted further to the high energy regime as 
compared with those seen in Fig. 3 or other panels in Fig. 4. Nevertheless, we still observe a 
clear window of almost zero conductance, but now with its location also significantly shifted 
as compared with the profile of Xfji^k). For the case of |?7 max | = 0.2w in Fig. 4(c), it is 
somewhat similar to the dotted line in Fig. 3(b), consistent with the fact that they have 
the same roughness strength. However, because here the location of the peak of Xfj{2k) 
is close to the threshold k value, the zero conductance window is also near this threshold: 
the conductance curve rises when k exceeds the threshold and then it quickly drops to zero 
again. For the case of |?7 max | = O.lio, its zero conductance window shown in Fig. 4(b) is 
narrower than those in Fig. 4(c) and Fig. 4(d), consistent with our intuition. Somewhat 
surprising is the case shown in Fig. 4(a), where the roughness strength is weak and the 
energy threshold for transmission is almost unaffected. But still, a narrow window for very 
small conductance is clearly seen in Fig. 4(a). This result is unexpected, because if one 
applies Eq. ( |T2|) directly, one would predict that no such conductance window should occur 
for l^maxl = 0.01. Further, the conductance window in Fig. 4(a) is much left-shifted as 
compared with the profile of Xni^k) (inset of Fig. 4(a)). Similar results are obtained in 
other realizations of the surface roughness function t](x) that have similar profile of the 
structure factor 22|]. Such a remarkable deviation from the theory, we believe, is due to 
a breakdown of the Born approximation in deriving Eq. ( |T2i) . Indeed, the conductance 
window for |r/ max | = 0.01 is located in a regime of very low scattering energy, and is hence 
not describable by a theory based on the Born approximation. Certainly, it should be of 
considerable interest to experimentally study the conductance windows in these cases of 
weak surface roughness. 

To further confirm that the conductance windows observed here are due to the colored 
surface disorder, we note that if we consider a surface function as that shown in the inset of 
Fig. 3(a), then all the conductance windows shown here indeed disappear and the results 
will become similar to what is shown in Fig. 3(a). 
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FIG. 5: Conductance properties of rough bended waveguides, with a straight upper boundary 
P(x) = 1 and a bended lower boundary given by Q(x) = 4a(L/2 — y) 2 /L 2 + fj(x). The waveguide 
curvature parameter a is given by a = in (a), a = 0.5 in (b), and a = 1.0 in (c). The inset in (a) 
shows the surface structure factor used in all the calculations, and the insets in (b) and (c) show 
the associated waveguide geometry. 

B. Bended Rough Waveguides 

In Fig. 5 we examine the conductance properties of a bended rough waveguide (Fig. 5(b) 
and Fig. 5(c)) as compared with those of a straight rough waveguide (Fig. 5(a)). In all 
the three cases shown in Fig. 5, the upper boundary is given by P(x) = 1; and the lower 
boundary is a parabolic curve plus random fluctuations, i.e., Q(x) = 4a(x — L/2) 2 / L 2 + fj(x); 
with a = in Fig. 5(a), a = 0.5 in Fig. 5(b), and a = 1.0 in Fig. 5(c). As to the structure 
factor of fj(x), it is assumed to be of a double-window form as shown in the inset of Fig. 5(a), 
with |r7 ma x| the same as in Fig. 3. In the case of a straight rough waveguide, this double- 
window structure factor creates an analogous double-window structure in the conductance 
curve (Fig. 5(a)), with its location matching the profile of the structure factor. Interestingly, 
as we introduce a curvature in the lower boundary in Fig. 5(b), the double- window structure 
survives but shifts considerably toward higher k values. In Fig. 5(c), the curvature of the 
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FIG. 6: (a) Conductance of a rough waveguide with symmetric surfaces, modeled by an upper 
boundary P{x) = 1 + fj and a lower boundary Q(x) = fj(x). Inset shows the surface structure factor 
used in all the calculations here, (b) Conductance of a rough waveguide with two uncorrelated 
surfaces, modeled by an upper boundary P(x) = 1 + fj'(x) and a lower boundary of Q(x) = fj(x). 
Inset shows the associated waveguide geometry, (c) Conductance of a rough waveguide with two 
antisymmetric surfaces, modeled by an upper boundary P{x) = 1 — r)(x) and a lower boundary 
Q(x) = fj(x). Inset shows the associated waveguide geometry. 

rough waveguide further increases, the transmission threshold value of k also increases (as 
expected) , and the fingerprints of the double- window structure factor can still be seen in the 
conductance curve. We have also checked if we create three windows in the structure factor, 
then three windows in the conductance curves can be induced as well, with their locations 
controllable by tuning the curvature of the bended waveguide. 

Finally we consider waveguides with both upper and lower boundaries being rough. In- 
terestingly, in this case a more sophisticated theory [14{| shows that the scattering can be 
regarded as the scattering in a smooth waveguide plus an additional effective potential. The 
theoretical electron mean free path, calculated using a Green function averaged over different 
surfaces, is shown to be contributed by different terms, due to different mechanisms called 
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amplitude scattering, gradient scattering, and square gradient scattering [1J|. The impor- 
tance of these terms depends on whether the upper and lower boundaries are symmetric, 
uncorrelated or anti-symmetric. Motivated by this interesting prediction, here we show in 
Fig. 6 three computational results, for symmetric (Fig. 6(a)), uncorrelated (Fig. 6(b)), and 
anti-symmetric boundaries (Fig. 6(c)), all with the same roughness strength as in Fig. 3. 

For the symmetric case, the effective waveguide width is not affected by the roughness. 
By contrast, for the antisymmetric case of Fig. 3(c), the effective waveguide width is most 
affected. These two simple observations explain why the threshold k value for transmission 
is the smallest in Fig. 6(a) and the largest in Fig. 6(c). Even more noteworthy is how 
the structure factor of surface roughness generates a conductance window for these three 
cases. In particular, the window of the conductance drop in the symmetric case (Fig. 6(a)) 
is narrower than that seen in Fig. 6(b) and Fig. 6(c). Moreover, the conductance window 
in the anti-symmetric case is the widest one, and is much shifted towards high k values as 
compared with the structure factor. This large mismatch between the conductance window 
and the peak location of Xrj(2h) hence reflects clearly an effect from the correlation between 
the two rough boundaries. Though our results cannot be easily explained by the theoretical 
result of Eq. (TT2l , they are consistent with the theoretical prediction in Ref . [14J that among 
the three cases of symmetric, uncorrelated and anti-symmetric waveguides, the electron mean 
free path in anti-symmetric waveguides should be the shortest. 



IV. DISCUSSION AND CONCLUSION 



In this computational study we have focused on how the structure of surface roughness 
impacts the conductance properties of electrons propagating in a quantum wire modeled 
by a 2D waveguide. Our conductance results are directly computed from a reaction matrix 
approach. An early theoretical result is hence confirmed by detailed behavior of the con- 
ductance, a quantity that should be measurable in experiments. In addition, our results for 
symmetric, uncorrelated, and anti-symmetric rough waveguides are consistent with a very 



recent theory 



Unlike in the bulk case, for quantum wires of limited length the sensitive dependence of 
the localization length upon the structure factor of surface roughness can be easily manifested 
in conductance properties. Our direct scattering calculations show that this is true, even 
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for those interesting cases that are beyond the domain of today's theory or have not been 
treated theoretically. We conclude that conductance properties are easily controllable by 
engineering the surface roughness of quantum wires. 

Though we have focused on the transport behavior of electrons, we believe that our 
methodology might be also useful for studies of other types of wave propagation in disordered 
systems. In particular, there is now a keen interest in understanding phonon transport in 
rough quantum wires. Recent computational work and experiment work [24]] showed the 
importance of surface disorder in the heat transport of thin silicon nanowires with a radius of 
w = 22 nm. It was also demonstrated experimentally that surface roughness can be used to 
dramatically suppress heat conductivity 6| and hence enhance thermoelectric efficiency for 
thin silicon nanowires with a radius about w = 50 nm. Our computational tools, together 
with the guidance from the theory 

aaaaa 

, might help answer some important 
questions regarding to phonon transport in rough nanowires. Indeed, we conjecture that 
it should be possible to design some colored surface to create conductance windows for 
phonons, but not for electrons. If this is indeed realized, then electron conductance is 
not much affected and phonon conductance will be greatly reduced. This will be of vast 
importance to thermoelectric applications. 

Finally, we note that spin accumulation in quantum waveguides with rough boundaries 
was recently studied in Ref . [3] . It should be interesting to see how colored surface disorder 
might have some useful impact on spin accumulation effects or spin transport. 
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